Kirchhoff 's Loop Law and the maximum entropy production principle 
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In contrast to the standard derivation of Kirchhoff 's loop law, which invokes electric potential, we 
show, for the linear planar electric network in a stationary state at the fixed temperature, that loop 
law can be derived from the maximum entropy production principle. This means that the currents 
in network branches are distributed in such a way as to achieve the state of maximum entropy 
production. 
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I. INTRODUCTION 

Kirchhoff's laws Jj are the standard part of general 
physics courses 0, U ■ In electrical engineering they 
are the starting point for the analysis of stationary pro- 
cesses in electric networks 0- In the stationary state, 
due to the principle of charge conservation, Kirchhoff's 
current law (current law) is valid. It states that in each 
node of the network the sum of ingoing currents equals 
the sum of outgoing currents. Kirchhoff' loop law (loop 
law) is based on the assumption that electric potential is 
a well defined quantity in any point of the electric net- 
work. Then one can apply the principle of energy con- 
servation to a macroscopic small amount of the charge 
circulating around the loop, i.e. the energy obtained on 
the sources should be equal to the dissipated energy. This 
statement is equivalent to the loop law, which states that 
the algebraic sum of electromotive forces (EMF-s) of the 
sources is equal to the sum of voltages (potential differ- 
ences) in the loop. In this paper we show that the loop 
law can be derived for a linear planar network using the 
maximum entropy production principle |^]. This means 
that stationary state currents distribute themselves in 
the branches in such a way as to maximize the entropy 
production in the network. 



II. MESH CURRENTS AND LOOP LAW 



We consider a planar network (see Fig. 1). Assuming 
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FIG. 1: A linear planar electric network 



that network parameters, EMF-s and resistances, are 
fixed, one can find all currents applying current and loop 
law to nodes and loops. However, due to the current law, 
the currents in branches are not independent quantities. 

Kirchhoff's laws give no prescription on how to find a 
set of independent currents for a given electric network. 
In the case of the planar network this problem has been 
solved by electrical engineers jsl by introducing the con- 
cept of mesh currents. Let us first define, within the 
network, the simple loop (mesh) as the one having no 
loop within it (loops 1,2 and 3 in FigCJ. We associate 
a mesh current with each mesh (Ji, J2 and J3 in Fig. 
A current in a branch, common to two neighbor- 
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FIG. 2: Mesh currents and Kirchhoff 's current law 



ing meshes, is an algebraic sum of corresponding mesh 
currents (see Fig|21). The current in the outer branch, 
the branch which belongs to one mesh only, is equal to 
the mesh current. Evidently the mesh currents incorpo- 
rate the current law (see Figl^J. It is easy to prove, by 
means of the mathematical induction, that the number of 
contour currents is equal to the number of independent 
currents in the network. The mesh currents are indepen- 
dent parameters determining the stationary state of the 
electric network as the thermodynamic system. In order 
to make the analysis of the network in terms of the mesh 
currents as simple as possible we introduce equivalent 
EMF-s and resistances. The equivalent EMF is equal to 
the algebraic sum of the EMF-s in a certain branch and 
the equivalent resistance is the sum of the resistances in 
that branch. We enumerate the meshes and correspond- 
ing mesh currents by single index notation, while the 
double index notations is used for the equivalent EMF-s 
and for equivalent resistances (Fig^. Different indices in 
the double index notation appear when the single branch 
is shared between two meshes. 

Applying the loop law for each mesh loop in terms of 
mesh currents we obtain the system of linear equations, 
the number of equations being equal to the number of 
mesh currents, 



(1) 



Note that the system of equations incorporates both 
Kirchhoff's laws (due to the definition of the mesh cur- 
rents). Here the left hand side of the equation is the 
algebraic sum of EMF-s within i ~ th mesh loop. The 
EMF is positive if mesh current comes out of its positive 
pole. Otherwise it is negative. Evidently it holds true 
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^ij — R'ji- 



(3) 



We stress that the system of equations ^ assumes that 
the principle of energy conservation holds for each mesh 
separately (loop law). In the following we show that this 
system of equations ^ can be alternatively derived ap- 
plying the maximum entropy production principle under 
the condition that the principle of energy conservation is 
valid for the whole network as the thermodynamic sys- 
tem. 



III. 



CONSERVATION OF THE ENERGY 



If the system at the fixed absolute temperature T re- 
leases heat per unit of time, dQ/dt, the corresponding 
entropy production is 



djS 
~dt 



1 dQ 
T~dt' 
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where index i, as it is introduced in the theory of 
non-equilibrium processes 0, emphasizes that only ir- 
reversible entropy change (rise) is considered. 

The entropy production is the function of the ther- 
modynamic state of a system. The stationary state of 
resistor at the fixed temperature is defined by the mag- 
nitude and the sense of the flow of the electric current. 
Due to the isotropy of the resistor the entropy produc- 
tion is invariant with respect to the flow direction of the 
current. 



(5) 



Assuming that the stationary state of the resistor is close 
to its equilibrium state Taylor expansion of the entropy 
production gives, 



^ = a + bJ + c.P + d.P + 0{J^). 
dt 



(6) 



A free term vanish since it corresponds to the entropy 
production in the equilibrium state. Due to the relation 
(O coefficients of all odd powers vanish while c > due to 
dSi/dt > 0. Putting c — R/T, where R is the resistance 
we get the well known expression for the dissipated elec- 
tric power dQ/dt = RJ^. Evidently c > implies R> 0. 

A linear planar network can be separated into sources 
of the forces (EMF-s) and the passive part (resistors). In 
the stationary state there is no change in internal energy 
of the resistors while they convert energy from the sources 
into the heat given off to the surroundings. The rate of 
the energy conversion is given by. 



dt 



(7) 
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Energy released by all sources (EMF-s) per unit of time The solution of the system of n equations 



IS 



dW 
~dt' 



'^E^^Ji + - '^Eij{Ji- Jj) (8) 



The factor i, in these equations, appears due to the dou- 
ble counting in the case of internal branches. 

Bearing in mind that resistors do not change their in- 
ternal energy we can write, according to the first law of 
thermodynamics. 



dQ dW 
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dt dt 

i.e. taking into account and it holds, 
EiiJi + 2 X! - Jj) ^ 

i ij 



(9) 



(10) 



IV. ENTROPY PRODUCTION AND ITS 
EXTREMUM IN THE ELECTRIC NETWORK 

We argue that Kirchhoff's loop law follows from the 
principle of the maximum overall entropy production in 
the network, assuming that the energy conservation H10() 
is satisfied. If there are n meshes, and n associated mesh 
currents, we have a conditional extremum problem in the 
n-dimensional linear space. 

Assuming that all resistors in the network are at the 
same temperature T the maximum of entropy production 
occurs at the same point of the n— dimensional space of 
currents as the maximum of generated heat (see eq. 

Standard procedure ^ of solving the conditional ex- 
tremum introduces Lagrange's multipliers. In this case 
one seeks the extremum of the function 



(11) 



where due to the condition (|10|l 

3 id 
i i,j 

The function F now reads 



0. 



(12) 
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that satisfies the condition (|12|l is represented by the 
point in the rt-dimensional liner space { Ji} in which the 
function dtS/dt exhibits extremum. The value of La- 
grange's multiplier A is uniquely determined by the sys- 
tem of equations (|14|l and eq. (|12|l . 

The standard, rather tedious, procedure of determin- 
ing A can be avoided in the following way. Let us multiply 
eq. H14|l with Ji and sum over all mesh currents, i.e. 



2(1 -A) 



'^RiiJf + '^Rij{Ji - Jj)Ji 

i ij 

+ A^^yJ, =0. (15) 



Due to the symmetry relations Q and (j3Jl we have, 

Y,E,,J, = -Y,E,3J3, (16) 

ij ij 

and 

ij ij 

By making use of eqs. (|16|) and (|17|l the ea. H15|l becomes. 



2(1 -A) 



+A 



^ + - ^ Rij{Ji - Jjf 



Y,E,,J, + ]-Y,E^,{J^- Jj) =0. (18) 



Substituting the term linear in mesh currents in eq. H18|l 
by the term quadratic in mesh currents according to eq. 
(inji we get 



(2 -A) 



'^Rujf + - ^i?y (Ji - Jjf 



0, (19) 



i.e. A = 2. One easily finds that the system of equations 
(|14|l for A = 2 is just the system of linear equations J^l, 
which expresses the loop law. 

Let us add that F function for A = 2 is the n- 
dimensional paraboloid put upside-down, i.e. its extreme 
is maximum. In this way the equivalence between the 
maximum entropy production principle and Kirchhoff's 
loop law is established. In other words we can say that 
currents in a linear planar network distribute themselves 
so as to achieve the state of maximum entropy produc- 
tion. 
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V. DISCUSSION 

It was Ehrenfest (Enzykl. Math. Wissensch, IV, 2(11) 
fasc.6, p82, note23, 1912.) who first asked the question 
if there exists sorae function which, like the entropy in 
the equilibrium state of an isolated system, achieves its 
extreme value in the stationary non-equilibrium state. 

The minimum entropy production theorem, attributed 
to Prigogine , identifies the entropy production due to 
irreversible processes in the system, as a physical func- 
tion which becomes extremal in a stationary state. We 
agree with Jaynes opinion that Prigogine's theorem is 
in fact the definition of the very special stationary state 
with zero induced flux^known as the static head state in 
electrical engineering Q and in bioenergetics I* is 

the quasi-equilibrium state that can be established close 
to equilibrium in the case of linear relationship between 
forces and fluxes, when induced force is left free to seek its 
maximal possible value. Energy conversion then stops at 
the level of driving force and corresponding flux, because 
induced flux vanishes and cannot be used to perform any 
work. One cannot obtain phenomenological equations, 
such as those found in references 0,^3 '^^ ^1- from 
Prigogine's theorem, because that theorem docs not have 
power of a principle. On the other hand phenomenolog- 
ical equations, such as eq. JQ), can be derived from the 
maximum entropy production principle. 

For the one-loop network the maximum entropy pro- 
duction principle holds as well. We leave for a reader to 
repeat our calculations from equations Ullfl to H19|l in the 
case of the one-loop network. The claim that entropy pro- 
duction is minimal in the stationary state of the one-loop 
network |llj | is wrong on two accounts. Firstly, authors 
[Tll | have presented potential drops as thermodynamical 
forces, i.e. they make no distinction between sources 
and loads and introduce artifical " forces" and " currents" 
in the entropy production expression. In the station- 
ary state the charge conservation law requires that only 
one current fiows in the one-loop network, which must 
be driven by single thermodynamic force (an algebraic 
sum of all EMF-s in the loop). Secondly, the statement 
about entropy production minimization in ref. |llj | im- 
plies that minimum entropy production theorem |3 holds 
for a such network. This is not possible, because a min- 
imum of two thermodynamic forces and corresponding 
currents and the condition of vanishing induced current 
(the static head non-equilibrium condition) must exist 
for the minimum entropy production theorem to hold. 

An arbitrary network with no sources of electromo- 
tive forces (EMF) was first considered with respect to 
the entropy production by Kirchhoff |p , Maxwell ^3 i 
and Jeans (see §357. in reference 13]). Kirchhoff pos- 
tulated variational theorem 9J, which states that in an 
arbitrary volume with fixed surface potential the currents 
distribute themselves so as to achieve the state of min- 
imum entropy production. Since this example contains 
no sources of EMF, it represents a continuous extension 
of discrete network with no sources which was considered 



by Maxwell j3| and Jeans (see §357. in reference |T3|). 
These authors came to the same conclusion that min- 
imum entropy production governs current distribution, 
because they did not consider the energy conservation 
principle. 

When EMF-s and the energy conservation principle 
are taken into account then the principle of the maxi- 
mum overall entropy production in the network leads to 
Kirchhoffs loop law and accordingly to observed station- 
ary current distribution as demonstrated in this paper. 
Let us note that Jeans considered an arbitrary, not only 
planar, electric network with EMF-s too. The maximum 
entropy production principle is implicitly present in his 
theorem (see §358. in reference 13] ). This means that 
the current distribution in non-planar networks is also 
determined by the maximum entropy production prin- 
ciple. Jeans did not define independent currents in an 
arbitrary network, as we have done by introducing mesh 
currents in full agreement with charge conservation law. 
Therefore, it is relatively hard to follow his derivation. 

Our macroscopic approach is very similar to that of 
Onsager who considered the related problem of heat con- 
duction in an anisotropic crystal ■ We shall show in a 
separate paper (in preparation) that an equivalence ex- 
ists between the principle of least dissipation of energy 
[T^ Hsj l and the maximum entropy production principle. 
Onsager was aware of the fact that the stationary state 
of heat conduction is in fact the state of maximum en- 
tropy production (see the sixth section in the reference 
[Tjl)- The system of equations that describes station- 
ary processes can be inferred from the maximum entropy 
production principle both in the case of heat conduction 
and in the case of linear planar electric network. The 
advantage of the maximum entropy production princi- 
ple is that it offers a better physical insight through the 
explicit use of the energy conservation law in the con- 
ditional extremum problem. Rate of work done must be 
equal to heat produced in a stationary state. The dissipa- 
tion function originally introduced by Lord Rayleigh [lq | 
and also used by Onsager [l^, obscures the distinction 
between energy sources and passive elements in the sys- 
tem, between the rate of work done and heat produced. 

Beside macroscopic formulation of the maximum en- 
tropy production principle there exists the microscopic 
formulation too. It is due to Kohler who demon- 
strated, starting from the Boltzmann transport equa- 
tion, that fiuxes in the stationary state of gases distribute 
themselves in order to achieve the state of maximum en- 
tropy production. Ziman |18| extended his work to free 
electron system in solids and demonstrated how this prin- 
ciple can be used in obtaining more accurate solutions 
of the Boltzmann transport equation. Recently, Dewar 
has shown, applying Jaynes information theory for- 
malism of statistical mechanics ^9|, |20j , that stationary 
states are characterized by maximum entropy produc- 
tion. 

To conclude, it is our belief that the principle of 
maximum entropy production is valid in linear non- 
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equilibrium thermodynamics, at least for stationary pro- 
cesses. In other words we are convinced that it is just 
the principle Ehrenfest was looking for. We stressed in 
this paper that in each situation (linear planar electric 
network in this paper) this principle must be applied as 
the conditional extrcmum problem so that energy con- 
servation law holds. Two of us have shown recently that 
biochemical cycle kinetics close to equilibrium state can 
be described by an analogue electrical circuit [U, and 
that modeling of photosynthesis ^2 can also be done by 
using the maximum entropy production principle. These 
and other older 0| and more recent 0, |2j| results using 



this principle show that the present derivation of Kirch- 
hoff's loop law in this paper is likely to be only the be- 
ginning of the widespread use of this principle in many 
different scientific disciplines. 
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